
LOAN COPY: RETURN TO 
AFWL (WLIL-2) 

KIRTLAND AFB, N MEX 

ANOMALOUS DIFFUSION IN A PLASMA 
FORMED FROM THE EXHAUST BEAM OF 
AN ELECTRON-BOMBARDMENT 
ION THRUSTER 

by Allan J. Cohen 

Lewis Research Center 
Cleveland, Ohio 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • 

i 



WASHINGTON, D. C. 


AUGUST 1968 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 



0131254 

NASA TN D-4758 


ANOMALOUS DIFFUSION IN A PLASMA FORMED FROM THE EXHAUST 
BEAM OF AN ELECTRON-BOMBARDMENT ION THRUSTER 

By Allan J. Cohen 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 - CFSTI price $3.00 



ABSTRACT 


A 10-cm-diam. cylindrical anode with a grounded coaxial center rod was immersed 
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in the low density (n g ~ 2x10 cm ) exhaust beam of an ion thruster. A longitudinal mag- 
netic field was applied. A theoretical and experimental study was made of the electron 
diffusion process in the beam plasma passing through the electrically floating anode. 
Anomalous diffusion and low frequency noise were found. The mechanism of anomalous 
diffusion was considered similar to that reported by Kadomtsev and Nedospasov. Bohm’s 
diffusion equation was found applicable and a constant of 1/8 (twice Bohm's value) was 
deduced experimentally. Rotation in the plasma was found by means of a correlation 
technique. The measured frequency of rotation ranged from 100 to 150 kHz. 
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SUMMARY 

A 10-centimeter-diameter, cylindrical anode with a grounded coaxial center rod was 
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immersed in the low density (n g ~ 2X10 cm" ) exhaust beam of an ion thruster. A longi- 
tudinal magnetic field ranging up to 1 10x10" ^ tesla was applied. A theoretical and ex- 
perimental study was made of the electron diffusion process in the beam plasma passing 
through the electrically floating anode. Anomalous diffusion in the radial direction and 
low frequency noise were found. The mechanism of anomalous diffusion was considered 
similar to that reported by Kadomtsev and Nedospasov. Bohm's diffusion equation was 
found applicable and a constant (twice Bohm's value) was deduced experimentally. Rota- 
tion in the plasma was found by means of a correlation technique. The measured fre- 
quency of rotation ranging from 100 to 150 kilohertz was predictable from theory. 


INTRODUCTION 

A weak magnetic field is used in the discharge chamber of the electron- bombardment 
ion thruster (ref. 1) to lower the value of energy required to produce a beam ion and 
correspondingly increase efficiency. The magnetic field reduces the radial transport 
rate of the electrons across the discharge chamber of the ion thruster and as a result 
enhances performance. It has been observed that the electron diffusion process in the 
ion thruster falls into two regions, stable and unstable (refs. 2 and 3), and that the cor- 
responding diffusion coefficient behaves classically and nonclassically, respectively. It 
also has been observed that the performance of the ion thruster is significantly enhanced 
by the magnetic field only when the thruster is operating in the stable region (ref. 2). 
According to reference 3 a critical value of the magnetic field separating stable and un- 
stable regions can be calculated (following the Kadomtsev and Nedospasov treatment, 
ref. 4) for the discharge plasma of the ion thruster. It is important to note that the dis- 



charge plasma not only is subjected to an axial magnetic field in the cylindrical ion- 
chamber configuration but also contains an applied radial electric field and a longitudinal 
electric field which develops from potential penetration from end wall plasma sheaths. 

One of the first successful plasma diffusion theories for determining the boundary be- 
tween stable and unstable, or anomalous regions of a plasma, was the Kadomtsev and 
Nedospasov theory (ref. 4). In this theory rotating perturbations of number density and 
potential are introduced in a Fourier expansion and values for the critical magnetic field 
and rotational frequency are obtained. The theory was applied to the high pressure 
(~ 1 mm) positive column in a longitudinal magnetic field. An applied axial electric field 
was also present. Guest, Simon, and Hoh (refs. 5 and 6) further expanded the theory to 
less dense plasmas and to plasmas with properties other than those of the positive column. 
Hoh's physical explanation of this type of anomalous diffusion (the Kadomtsev and 
Nedospasov type) was to suggest that a density concentration of electrons revolves in the 
plasma moving through the neutral background and that an ion concentration follows 
closely behind, but slightly separated, thus giving rise to a rotating azimuthal potential 
gradient. As a result a Hall drift exists, when the magnetic field strength is greater than 
the critical value, enhancing motion across the magnetic field and causing an enhanced 
loss from the "classical" containment characteristic of the magnetic field. Since then 
several investigations of the diffusion process in a magnetized plasma have revealed ro- 
tation (see, for example, ref. 7) and many investigations have shown increased or anoma- 
lous diffusion coefficients (ref. 8). It is the belief of some investigators (ref. 8) that any 
partially ionized plasma submerged in a magnetic field and having an electric field either 
applied or developed is susceptible to anomalous diffusion beyond some critical magnetic 
field or beyond a critical value of a related parameter. 

Classically, the collisional diffusion coefficient in a magnetic field should behave as 
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B" (ref. 9). Early in plasma physics experimentation, Bohm (ref. 10) made a measure- 
ment of the diffusion coefficient in a magnetic field and found that the measured value dif- 
fered significantly from the classical value. He proposed a new empirical relation for 
the diffusion coefficient which contained a B” 1 variation with magnetic field. He believed 
this to be the proper variation of the diffusion coefficient for a plasma undergoing anoma- 
lous or, as he termed it, drain diffusion. Since that time numerous experiments (ref. 8) 
and theories (refs. 11 to 13) have been put forward to clarify the diffusion process in a 
magnetized plasma. 

Because of the inherent difficulties in plasma physics experiments aimed at deter- 
mining transport rates, the proper value of the diffusion coefficient in a magnetized 
plasma is not clear, nor is the physical explanation of its unstable behavior. It is the 
author's belief that anomalous diffusion in a magnetized, partially ionized plasma with 
electric fields can be understood in terms of the mechanism used by Kadomtsev and 
Nedospasov to predict stable and unstable regimes of particle transport. Throughout this 
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report, reference to the Kadomtsev and Nedospasov theory should be taken in the broad 
sense to include all those theories of anomalous diffusion in partially ionized gases in 
which there are rotational number density and potential perturbations and in which the 
pertinent mechanism of transport of the perturbations is the Hall effect. This mechanism 
can produce the Bohm diffusion that appears to predominate in the unstable mode of the 
magnetized plasma (see refs. 13 to 16). 

In the present report the electron diffusion process in a magnetic field is studied by 
floating a metal anode in the exhaust beam of an ion thruster. This configuration is 
studied instead of the discharge chamber process of the ion thruster because of its rela- 
tive simplicity. Although an exact correspondence is not expected, certain inferences 
are drawn with respect to the diffusion process of the ion thruster (and also the diffusion 
process in general) based on similarities. Because the metal anode is floated, a radial 
electric field develops and electrons drift radially in the combination electric and mag- 
netic field. 

The plasma is probed for electric potential and number density gradients to deter- 
mine whether the diffusion process is classical or anomalous in the range of magnetic 
field strengths studied. The Kadomtsev and Nedospasov theory is compared to the re- 
sults obtained. The plasma is also investigated for various values of electron current 
flow to the anode as well as electric field strength. Bohm's equation for the mobility 
coefficient with its characteristic B" * magnetic field dependence is then compared with 
these data. Finally, the plasma is studied for rotation in the unstable region by a corre- 
lation technique. As a result of this work and previous studies by other authors, conclu- 
sions are drawn about the diffusion process in a magnetic field and about the behavior of 
the ion thruster. 


THEORY 

Electron Diffusion Process 


Sketch (a) shows the experimental arrangement; a schematic is given in figure 1(a). 

A cylindrical coordinate system with orientation as shown in sketch (a) is used throughout 


Plasma 

beam 
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the report. An electrically floating cylindrical anode and coaxial grounded rod are set in 
the exhaust beam of an ion thruster. A uniform axial magnetic field is applied. The ion 
beam is neutralized by secondary electrons from the tank walls and other grounded equip- 
ment, forming an approximately neutral plasma for the experiment. The ion number 
density of the plasma n at any point can be obtained from the relation 


J 


B 


= neU 


zi 


( 1 ) 


where J 0 is the beam current density measured in this study by a molybdenum probe 
and U . is the exhaust velocity (a constant) of the ion thruster. (Symbols are defined in 
the appendix. ) The plasma beam not only has an axial velocity but also has a radial velo- 
city component, that is, it spreads after leaving the ion thruster. As a result the anode 
is struck by the ion beam and, since it is electrically floating, it builds up a positive po- 
tential with respect to the grounded rod. The electrons in the plasma are pulled to the 
anode by this electric field and effectively neutralize it. The axial magnetic field impedes 
the electron flow to the anode. 

The radial current density of electrons in the experiment can be deduced from the 
equation 


3n sir 


( 2 ) 


where the first term is the current due to a number density gradient and the second term 
is the current due to a potential gradient. This equation with classical values of D 
and ju is used to determine if the radial electron transport rate is classical or anoma- 
lous. The equation, with 0-dependent terms, is given in reference 5. If the radial num- 
ber density gradient is neglected (experimental justification will be discussed later), the 
equation for electron current density across the system reduces to the following Ohm’s 
law expression: 


J re =n e^ie 


av 

dr 


(3) 


From Einstein's relation (ref. 8), the result is 


kT 


D le " “xe 


(4) 


where D and n are the electron diffusion and mobility coefficients, respectively, 
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and kT e /e is the electron temperature in electron volts. Equation (4) shows that the mo- 
bility and diffusion coefficient are directly proportional for a constant electron tempera- 
ture. 

Classically, from reference 5, the electron mobility coefficient (for sufficiently 
large values of <^ e T e ) is 


M 


ie 


1 

_ 1 / m e 

B(x> T 
e e 

B 2 \ eT e 


(5) 


where = eB/m e , r g = l/N 0 (crv) , N q is the neutral density (assumed to be much 
greater than the plasma number density n), and (av) is the velocity times the collision 
cross section averaged over the electron velocity distribution. Values of u> e T e for the 
experiment are given in table I which also lists other important experimental constants. 

Bohm, however, has suggested (ref. 10) that, for anomalous or drain diffusion, the 
diffusion coefficient is 


D 


le 


kT e C B 


B 


( 6 ) 


where B is the magnetic field and the constant according to Bohm is 1/16. Com- 
paring equation (4) with equation (6) results in the corresponding (Bohm) mobility 




le 



(7) 


Equation (5) can be compared with equation (7) to show the difference between the classi- 
cal electron mobility and the Bohm mobility. 

In the anomalous diffusion regime, according to the Kadomtsev and Nedospasov dif- 
fusion theory, number density concentrations drift across magnetic field lines as a re- 
sult of the Hall effect. In reference 13 a derivation of Bohm's equation is presented 
based on this (Hall) effect. In this derivation it was necessary to assume that the derived 
expression for the Bohm coefficient was independent of magnetic field. Reference 7 added 
experimental support to the derivation of reference 13. 

In the present report, equation (3) is used empirically in the anomalous region to ob- 
tain a mobility coefficient that is equated to the Bohm mobility (eq. (7)) to arrive at a 
value for the Bohm constant Cg. 
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Number Density and Potential Gradient Behavior 


The longitudinal and radial behavior of the number density and potential to be ex- 
pected in the present experiment can be deduced from the continuity equation, the elec- 
tron equation of motion, and appropriate assumptions. The continuity equation for the 
ions in steady state is, in vector notation, 


V- (n.U.) = 0 


(8) 


and for the electrons is 


v.(„ e u e ) = o 


(9) 


where U = (l/n)J and U refers to the drift velocity. Equation (3) is the equation of 
motion for the electrons. For application to the present experiment, the following sim- 
plifying assumptions are considered appropriate: 

(1) Axial symmetry is assumed. 

(2) Plasma neutrality is assumed such that n = n e = m. 

(3) The quantity dV/dr is a function only of r. 

(4) U zi is a constant. 

(5) n is a function only of z. 

(6) U ri is a function only of r. 

(7) 1J is a function only of z. 

(8) No perturbations arise from influence of the ion beam outside the anode. 

Writing out the continuity equations in scalar form with the electron equation of motion 
(eq. (3)) substituted in them yields 


ie v dr 


(rn™) + ± 
\ e dr / dz 


( n e U ze> = 0 


( 10 ) 


- — ( rn i U ri ) + — ( n i u zi) = 0 

r 9r ^ ri 9z 1 21 


( 11 ) 


These partial differential equations can be easily separated using the previous assump- 
tions and the separating constants K' and K T ' to yield the following results: 



(12) 
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( 13 ) 


a_V _ K"r + ^1 
3r 2 Mie r 


n = n Q e 


(K'/U zl )z 


»n„ e 


-Cz 


(14) 


and 


U ze = 


K’ 


.Cz 


U zi + C 2 


n r 


(15) 


The constant C ^ is determined from the boundary condition that the electron current 
to the metal rod on the axis is zero, because there is no electron emission from the rod, 
and the plasma electrons are attracted to the anode. Thus, from equation (3), SV/3r is 

O 

zero at r = r j, the radius of the rod, so that = -K"r j/2jUj_ e * Equation (13) then be- 
comes 


3V 

9r 



(16) 


(Sheath effects are neglected in this study because, from table I, the Debye length is 
small compared with the dimensions of the experiment. ) Since the rod radius is small 
compared with experimental radii, the second term in equation (16) and hence also in 
equation (13) can be neglected. Integrating equation (13) then yields 


V = 


a l r 


+ b. 


(17) 


where a^ and b^ are constants. 

As a result of the preceding analysis, the behavior of the number density (eq. (14)) 
and the behavior of the electric potential (eq. (17)) have been deduced. As will be seen 
in later sections, the exponential in equation (14) was such that, for the experiment, the 
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behavior of the number density differed only slightly from a linear decline with in- 
creasing z. 


APPARATUS AND PROCEDURE 
Experimental Description 

The experimental setup is shown in figure 1. A collimated ion beam supplied from 
an ion thruster passed through an electrically floating 20- centimeter- long by 5-centimeter- 
radius cylindrical metal anode with an electrically grounded axially centered metal rod. 

The anode was electrically connected to a vacuum tube voltmeter of very high input im- 
pedance (see sketch (b)), which essentially measured the floating potential of the anode 


Anode 


Metal rod 


VTVM 


X 



(b) 


with respect to ground. Floating potential data were recorded for various magnetic field 
strengths. A microammeter was used in the metal rod leg to check current flow to this 
electrode. A small current related to ion impingement was found. 

The anode was set inside an axial magnetic field with a field variation of less than 
5 percent. Movable molybdenum probes, a description of which is given in the following 
section, were set at the entrance and exit of the anode and could be moved across the 
anode diameter. Two movable Langmuir probes were placed inside the anode and were 
free to move in the anode across a radius. The Langmuir probes were set at the center 
plane of the anode (z = constant) and were displaced by 120°. The motion of the probes 
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was obtained through an actuating mechanism driven by a low rpm motor. The probe po- 
sition was obtained from a potentiometer connected to an ohmmeter calibrated in linear 
distance. 

Typically the ion thruster accelerated to 2000 volts a 50-milliampere beam of mer- 
cury ions which were directed along the axis of the experiment. The beam was colli- 
mated and limited by a series of baffles before passing through the anode cylinder. 

About 2 milliamperes of beam current entered the anode, a significant part of which 
(about 50 percent) impinged on the anode wall. The ion thruster and experimental rig 
were placed in a vacuum chamber 1^ meters in diameter and 5 meters long which oper- 
ated at a pressure of lxlO" 6 to 2x10" 6 torr (1. 3xl0“ 4 to 2. 6xl0 -4 N/m^) when the thrus- 
ter was on. For calculating the classical electron mobility coefficient, it was found con- 
venient, however, to deliberately leak in air raising the pressure to 8x10" ^ torr 
(12. 4x10" N/m‘ ), for with this pressure it can be assumed that the neutral background 
consists mainly of nitrogen (see table I). Also, with this pressure the plasma potential 
fell between 10 and 75 volts, a convenient range for biasing the Langmuir probes. 

Finally a shielding screen was placed slightly beyond the exit of the anode to prevent 
plasma from contacting the outside surface of the anode cylinder (see fig. 1(a)). 


Molybdenum Probe Operation 


Sketch (c) shows the molybdenum button probe circuitry. The molybdenum button 
probe consists of a molybdenum disk about 3 millimeters in diameter facing into the ion 
beam. The probe is biased negatively by about 20 volts with respect to plasma potential 
to suppress electron collection. The impinging ion beam current is registered on an 
X-Y plotter by means of a voltage drop across a 10 000-ohm input resistor and plotted 
against probe position. The beam current density is obtained by dividing probe current 
by the area of the disk. The plasma number density can be obtained from equation (1) 
and the assumption of quasi- neutrality. Values of number density were obtained in this 
manner at the entrance and exit of the anode, and values of number density in the center 
plane were obtained by calculations based on equation (14). As was noted previously, 


Probe (— | 


Input <£ 
resistor *“> 


1 X-Y Plotter 
Bias | 
battery i 



Position 

voltage 
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equation (14) differed only slightly from linear in the present study. In addition to num- 
ber density data the molybdenum button probes can be used in determining anode ion im- 
pingement current, which is equal to electron current for an electrically floating anode. 
This is simply accomplished by integrating ion current density over the anode entrance 
area and exit area and then subtracting the ion current leaving the anode from the ion 
current entering the anode. The Langmuir probes, when operated in the mode shown in 
sketch (c), can also give number density data by taking into account the probe area sub- 
ject to ion impingement. 


Langmuir Probe Operation 

Sketch (d) shows the circuit of the Langmuir probe used in the experiment. 

An X-Y plotter was used in conjunction with the Langmuir probe to obtain electron or 
ion current to the probe as a function of bias potential. A 10 000-ohm input resistor was 
used. Typical Langmuir probe data for this experiment is indicated in sketch (e). In 
sketch (e) probe current Ip is plotted against bias potential V^, and the floating poten- 




(e) 


tial V, is identified as the bias potential corresponding to zero probe current. Typical 
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ranges of Ip and are 0 to 2X10 ampere and 0 to 50 volts, respectively. In this 
experiment the Langmuir probe was used to measure radial electric fields. To obtain 
this quantity, it was necessary to make measurements of plasma potential against radial 
position. From reference 17 plasma potential can be obtained from floating potential and 
electron temperature by the following relation, 

kT 

v pp - v t + g -f- <18 > 


10 




where V is plasma potential and g = 3 for mercury ions according to reference 17. 

The electron temperature can be obtained from the slope of a semilog plot of Langmuir 
probe data using the conventional probe analysis technique. In the experiment, plasma 
potential was measured at three equally separated radial positions. The inner and outer- 
most potential values were used with equation (17) to obtain the constants a^ and b^. 

The middle potential measurement is used as a simple check of the equation. By differ- 
entiation of equation (17), the radial electric field can be obtained. 

The aforementioned data (both molybdenum button and Langmuir probe data) were 
used in equation (2) to determine whether the diffusion process was classical or anoma- 
lous. Molybdenum button probe data yielded values of number density and also ion im- 
pingement current to the anode (equal to electron flow), while Langmuir probe data yielded 
measurements of radial electric field. An empirical value of then was obtained 

from equation (3) with J re , n, and 9V/3r being taken at the radius of the anode a, and 
in the center plane of the anode, where J re (a) was assumed to be an average value. This 
value of ju was compared to Bohm's empirical relation (eq. (7)) for different values of 
magnetic field. 


Correlation Technique 

According to the theory of anomalous diffusion considered in this report (the 
Kadomtsev and Nedospasov theory), number density perturbations as well as potential 
perturbations rotate in the anode cylinder. Both these types of perturbations passing a 
Langmuir probe will result in an alternating current signal. The signal picked up by the 
Langmuir probes when the plasma was experiencing anomalous diffusion was found to be 
a spectrum of many frequencies including a great deal of noise. To pick out the signal 
corresponding to the rotation (which was believed to be embedded in the noise), a corre- 
lation technique was used. Tuned circuits were coupled to two Langmuir probes separa- 
ted by 120° in the central plane of the anode. Sketch (f) below shows the circuitry in- 



(f) 
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volved. Both circuits were sharply tuned to the same frequency, the frequency of rota- 
tion, and allowed a small frequency band (a band width of about 2 percent of the tuned fre- 
quency) consisting of noise and rotation to pass. The signals from the two probes were 
then applied to a dual beam oscilloscope and one of them triggered both beams. The ef- 
fect of triggering was to allow many sweeps to fall on top of each other and yield a clear 
and well-defined trace for the signal controlling the triggering. The shape of the other 
trace depended on whether or not it had a definite phase relation with the controlling 
signal. If the second trace contained the rotational signal, then a clear second signal 
should be obtained on the oscilloscope (since it has a definite phase relation with the con- 
trolling signal) with 120° phase shift with respect to the triggering signal. If, however, 
the signals passing through the tuned circuits were just noise, then the second trace 
would appear to have a random relation with the first trace. The rotational signals were 
distinguished from the noise in this manner. 


RESULTS AND DISCUSSION 
Floating Potential and Plasma Noise 

As a result of electrically floating the anode in the exhaust beam of the ion thruster, 
a positive potential was obtained on the anode with respect to ground. A typical curve of 
these potentials as a function of magnetic field is shown in figure 2. These values of po- 
tential are below 100 volts. With a Langmuir probe placed in the plasma, the anode elec- 
trically floating, and a magnetic field applied, noise such as depicted in figure 3 is ob- 
served on an oscilloscope connected to the Langmuir probe. This noise, which is deduced 
to be made up predominantly of a low frequency, is not observed when the anode is 
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grounded (no electron transport) or when the magnetic field strength is less than 20x10 
tesla. These observations and their relation to the diffusion process will be discussed 
more fully in a later section. First, though, molybdenum probe data and Langmuir probe 
data will be analyzed and discussed. Anomalous diffusion, which was found to exist in the 
experiment, will be then discussed in relation to the aforementioned observations. 


Molybdenum Probe Data 

Current density profiles of the ion beam plasma at the entrance and exit of the anode, 
obtained with molybdenum button probes, are presented in figures 4 and 5. Also pre- 
sented in figure 6 is the beam plasma profile at the center plane of the anode obtained with 
a Langmuir probe operated as a molybdenum button probe. The data presented are for 
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various values of magnetic field strength and for a floating anode. A most interesting 
feature of these profiles is their departure from flatness at magnetic fields greater than 
zero. An explanation for the dip in the center (about r = 0) of the back profiles might be 
that the grounded center rod supporting structure shielded some of the beam from im- 
pinging on the molybdenum disk. 

The bell shape of the front profiles and the double peaks of the other profiles (ob- 
served if the traces are reflected through r = 0) which occur with magnetic fields greater 
than zero is not easy to explain since the ion beam should not be affected by these weak 
magnetic field strengths. It was further observed, however, that when the anode is elec- 
trically grounded the profiles are independent of magnetic field and thus are relatively 
flat. Note that by itself the radial electric field associated with the floating anode poten- 
tial should not be strong enough to affect the profile. Since the behavior of the plasma is 
stable when the anode is electrically grounded and since experimentally the profile is af- 
fected by whether or not the anode is floating or grounded, perhaps the changing shapes 
with magnetic field is linked to the anomalous mode of the plasma; the complete explana- 
tion however, has not been deduced. 

From entrance and exit molybdenum button probe measurements, values of ion cur- 
rent into and out of the anode are calculated by area integration and are presented in 
table n. Also presented are values of electron current to the anode (equal to the ion im- 
pingement current to the anode), values of C from equation (14), and values of number 
density at the center plane of the anode. The electron current to the anode is obtained 
simply by subtracting ion current out from ion current in. Electron number density at 
the entrance and exit of the anode is obtained (assuming n. ~ n g ) from equation (1) by 
using the exhaust velocity of the ion thruster and average values of ion current density. 
Values of number density at the center plane of the anode are obtained by using the afore- 
mentioned number densities and equation (14). 

/ -CZ 

For the experimental values of C in the quantity e , the number density shows 
an almost linear decline from entrance to exit of the anode. Impingement current to the 
anode (equal to electron current) and number density at the center of the anode are shown 
as a function of magnetic field in figure 7. 

Langmuir Probe Data 

Langmuir probe traces were taken on an X-Y plotter at three equally spaced radial 
positions in the center plane of the anode for three values of magnetic field strength. 

The radius of the anode is 5 centimeters and, starting from the center rod, Langmuir 
probe traces were taken at a 1. 25, 2. 50, and 3. 75 centimeter radii. A typical Langmuir 
probe trace of current to the probe as a function of bias voltage is shown in figure 8. 


I _ 
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This is a replot of a trace from an X-Y plotter in which there was noise variation of 
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±0. 05x10 ampere. From this and similar plots, values of floating potential and 

Maxwellian temperature were obtained. Figure 9 is a replot of figure 8 on semilog paper. 

From the slope of this plot, the electron temperature was determined to be 5. 1 electron 

volts. Electron temperatures for the other traces ranged from 3 to 6 electron volts. 

From values of floating potential and electron temperature, values of plasma potential 

were obtained using equation (18). Three data points were obtained for each magnetic 

o 

field strength. These are plotted in figure 10, and curves of the form a 1 r + bj (see 
eq. (17)) fit to the end points of each set of three are drawn. The center plasma potential 
points fall close to the theoretical curves in each case. Resulting values of a^, bp and 
the calculated values of electric gradient at the anode radius are given in table III. Study 
of the plasma potential plots of figure 10 and the floating anode potential curve of figure 2 
indicates that there is a relatively small voltage drop in the plasma compared to the large 
sheath drop at the electrically grounded center rod. 


Anomalous Diffusion and the Bohm Equation 

The classical value of the transverse electron mobility is calculated from equation (5) 
and plotted as a function of magnetic field in figure 11. The calculation is based on the 
average value of electron velocity, the average cross section, and a neutral background 
of nitrogen at the pressure of 8x10”° torr (12. 4xl0 _<a N/m 4 ) mentioned in the section 
APPARATUS AND PROCEDURE. 

Equation (2) is the classical (collisional) equation of electron current flow to the 
anode in the experiment. Using the number density gradient calculated from figure 6 in 
the center plane at the anode radius and the classical diffusion coefficient, it is found that 
the first term in equation (2) is orders of magnitude smaller than the experimentally 
determined value for J rg (a). For this reason, its contribution can be neglected. Equa- 
tion (3) was then used to determine whether the experimental diffusion process was clas- 
sical. Substituting into equation (3) average values of electron current density, center 
plane values of number density, and experimental values of electric field gradient yielded 
computed values for the electron mobility. Comparison of these values with the theoreti- 
cal values plotted in figure 11 showed that experimental values were much higher. For 
example, for a magnetic field of 55x10”^ tesla the calculated value of the mobility was 

20 square meters per volt per second whereas from figure 11 the theoretical value is 
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5X10” square meters per volt per second. Similar differences in the electron mobility 
coefficient existed for the other magnetic field values investigated. These results are 
orders of magnitude different and strongly suggest that anomalous diffusion is occurring 
in the range of magnetic fields in the experimental study. Further evidence of anomalous 
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diffusion can be obtained if equation (3) is solved for electric field using theoretical values 
of collisional mobility and experimental values of number density and current density. 
Resulting values of electric field are in the range of 10 4 and 10 5 volts per meter, which 
would indicate a floating anode potential of thousands of volts across the 5-centimeter 
gap. This is inconsistent with the floating anode potentials of figure 2. 

As previously stated, plasma noise was detected with a floating Langmuir probe when 
the anode was electrically floated and a magnetic field was applied. The frequency of the 
noise was determined to range from 0 to 200 kilohertz per second. The noise was unob- 
servable with the same instrumentation setup when the anode was electrically grounded or 

_4 

when the magnetic field dropped below about 20x10 tesla. It has been noted in many 
studies (see, for example, ref. 10) that low frequency plasma noise occurs in conjunction 
with anomalous diffusion and further that the onset of anomalous diffusion was also the 
onset of low frequency noise in the same frequency range observed here. This would in- 
dicate (1) that the onset of anomalous diffusion is about 20x10” 4 tesla in this study and 
(2) that anomalous diffusion did not exist for an electrically grounded anode. Both con- 
clusions are consistent with the Kadomtsev and Nedospasov theory. 

First anomalous diffusion occurred above 20X10” 4 tesla which means transition to 

anomalous diffusion in the experiment occurs in the "tens of gauss" range. In refer- 

2 2 

ence 4, the pressure involved was about 1 torr (1. 3x10 N/m ), and the critical magnetic 

p. 3 

field occurred in the "thousands of gauss" range; in reference 6, the pressure was 10 
12 

torr (1. 3x10” N/m ) and the onset occurred in the "hundreds of gauss" range. In the 
present experiment, the pressure is about 10” torr (1. 3x10” N/nr ) (see table I), and 
the onset occurring in the "tens of gauss" range is consistent with the general concept of 
lower critical fields with lower pressure. Finally, grounding the anode probably destroys 
the radial electric field (see eq. (18) with K = 0) and anomalous diffusion, which depends 
on electrons drifting in this electric field, therefore does not exist. 

Empirical values of electron mobility for three different values of magnetic field 
were obtained as before from equation (3) by substituting into it experimentally deter- 
mined values of electron current density, number density, and electric field gradient. 
These values are plotted in figure 12. The object is to see if Bohm's diffusion equation 
(eq. (7)) is applicable. A plot of Cg/B is drawn in figure 12 for Cg taken as 0. 12 
(which was found to be a good fit to the data). It can be concluded that the mobility be- 
haves as B”* (but with a coefficient Cg two times greater than Bohm's value). In this 
respect Bohm's equation is applicable to the data obtained herein. In reference 15 it was 
also found that the coefficient Cg was about 2 times greater than Bohm's value, al- 
though other values have been found by other investigators (e. g. , ref. 14). 
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Rotational Study 


Rotation of a number density concentration and potential gradient perturbation is 
intrinsic to the Kadomtsev and Nedospasov approach to instability in a plasma with a 
neutral background. Both conditions can induce an alternating current signal on a 
Langmuir probe placed in the plasma. A correlation technique described in a preceding 
section was used to detect rotation in the beam plasma. The plasma in the device studied 
in reference 7 was found to rotate with a velocity equal to one-fifth the velocity deter- 
mined by the radial electric field |(l/5)(E r /B)j . The plasma in the device studied in 
reference 18 rotated with a velocity equal to that determined by the radial electric field 
E r /B when the magnetic field was greater than the critical field. The rotational fre- 
quency for a rotational velocity of E r /B is 



for the m = 1 mode usually considered. From equation (17) this frequency is independ- 
ent of r. 

The correlation technique was used to determine whether rotation could be found in 
the frequency range of equation (19) (i. e. , 105 to 150 kHz/sec for the magnetic and elec- 
tric field range in the experiment). First, however, the tuned circuits were set at 
88 kilohertz per second, outside the rotational frequency range of the experiment, and the 
magnetic field was varied. A typical result is shown in figure 13. The lower signal in 
this figure as well as in all subsequent figures is the triggering signal. The lack of a 
clear upper trace in figure 13 indicates no significant correlation. The frequency was 
then adjusted to 125 kilohertz per second. Figure 14 resulted when the magnetic field 

4 

was set at 45x10 tesla. Correlation is apparent with the larger signal on top out of 
phase with the lower signal by about 125°. This is quite close to the 120° phase differ- 
ence in the Langmuir probe positions. The lower amplitude signal visible on top appears 
to be the result of another mode in the plasma rotation. Next, the tuned circuits are set 

_4 

at 140 kilohertz per second and correlation is found at a magnetic field of 97x10 tesla 
as shown in figure 15. The phase shift here is about 136°. Finally with the same set- 
tings the magnetic field leads are reversed and figure 16 resulted (a slight amount of 
returning was necessary). In this figure the upper wave is shifted forward, by about 
117°, instead of backward as in the previous figure, in just the manner expected for ro- 
tation under the influence of a reversed magnetic field. 

A plot of equation (19) is shown in figure 17 with frequency of rotation as a function 
of magnetic field represented. Values of electric field required in equation (19) were 
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taken from table IQ. Plotted also, as two experimental points, are the correlation data 
from figures 14 and 15. The points lie close to the curve consistent with equation (19) as 
the proper rotational frequency. 


Anomalous Diffusion with Respect to the Ion Thruster 

In reference 3 a calculation of the Kadomtsev and Nedospasov type was made to 
determine the boundary between stability and instability of the ionization chamber of an 
ion thruster. There was indication of anomalous diffusion in the ion thruster beyond a 
critical magnetic field in the form of noise measurements in the discharge chamber as 
well as in discharge chamber efficiency (ion loss per beam ion as a function of magnetic 
field). In the present experiment more direct plasma measurements were made from 
which anomalous diffusion was deduced. It is interesting to compare the two situations. 
Both cases deal with partially ionized plasmas. Both have cylindrical configurations with 
longitudinal magnetic fields, and both have associated radial electric fields in the plasmas. 
Rotation has been found in the present experiment but has not yet been definitely found in 
the ion thruster chamber. Both cases show an onset of low frequency noise and are 
anomalous in the noisy region (ref. 3). 

The similarities between the present experiment and the ion thruster discharge 
chamber suggest that the anomalous diffusion predicted in the ion thruster may result 
from the same type of mechanism as that obtained in this experiment. If true, this sug- 
gests that stabilizing magnetic field configurations (ref. 19) should be investigated as a 
possible means of reducing these loss mechanisms in the ion thruster discharge chamber. 


CONCLUDING REMARKS 

An electrically floating anode was set in a neutralized ion beam accelerated from an 
ion thruster. A theoretical study was made to determine whether the transport character- 
istics were of the classical or anomalous mode for the configuration. Floating anode po- 
tential data, molybdenum button and Langmuir probe data as well as correlation data. 
were taken. Furthermore, low frequency noise from a Langmuir probe was observed for 
various magnetic fields and for a floating anode condition. As a result of this study, 
anomalous diffusion was found, and its onset occurred at a low magnetic field strength. 

A relation between anomalous diffusion of the Kadomtsev and Nedospasov type and Bohm 
diffusion equation was thought to be linked by the Hall effect. Bohms equation was found 
to agree with experiment if the Bohm coefficient was increased by a factor of 2. Since 
the Kadomtsev and Nedospasov approach appeared applicable, a correlation study was 
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undertaken in order to find rotation. Rotation was found at the predicted values ranging 
from 100 to 150 kilohertz per second. The report supports many previous studies in 
predicting the onset of anomalous diffusion in a plasma with crossed magnetic and elec- 
tric fields. 

It was speculated that the electron diffusion in an ion thruster discharge may be the 
same type as that observed herein. If so, the loss might be alleviated with a stabilizing 
magnetic field configuration. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 20, 1968, 

120-26-02-10-22. 
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APPENDIX - SYMBOLS 


r| 

$ 


a 

anode radius 

U zi 

ion beam velocity 

a i' b i 

constants in eq. (17) 

V 

electric potential 

B 

magnetic field strength 

V 

velocity 

C 

K '/°zi 

X 

x-coordinate on oscilloscope 

C B 

Bohm's constant 

y 

y-coordinate on oscilloscope 

C D 

constant 

z 

z -coordinate 



6 

0-coordinate 

Ci,c 2 

constants 





x e 

Debye distance 

D ie 

electron transverse diffusion 

^ le 

electron transverse mobility 


coefficient 

E 

electric field 

a 

cross section 

e 

electron charge 

T e 

electron collision time 

I 

current 

w e 

electron cyclotron frequency 

J 

current density 

Subscripts: 

J B 

beam current density 

b 

bias 

K' , K” 

separation constants 

e 

electron 

l 

length 

f 

floating 

m e 

electron mass 

i 

ion 

N o 

neutral number density 

P 

probe 



PP 

plas ma 

n 

plasma number density 

r 

r- coordinate 

n o 

n(z = 0) 

t 

total 

r 

r- coordinate 

z 

z-coordinate 

i\ 

radius of centered rod 



l 


e 

0- coordinate 

T e 

electron temperature 



U 

drift velocity 
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TABLE I. - CONSTANTS IN THE EXPERIMENT 


-3 

Neutral number density, N rt , m 

17 

2x10 ' 

— ° 2 

Average cross section, a, m 

14x10" 18 

Electron collision time, t q , sec 

36X10"° 

Electron temperature, kT g /e, eV 

3 to 6 

n 

Debye distance, m 

0. 5x10" ^ 

o> e T e (B = 10xl0" 4 tesla) 
co e T e (B = 100x10” 4 tesla) 

6500 

65 000 


TABLE H. - COMPUTED RESULTS BASED ON 


MOLYBDENUM DATA 


Magnetic field, 
B, 

tesla 

Ion beam 
current, 
mA 

Electron 
current 
to anode, 
mA 

c ,a 

l/m 

Number density, 

n, 

m" 3 

In 

Out 

OxlO " 4 

1. 54 

0. 77 

0. 77 

2. 76 

1. 88X10 13 

b 28 

— 

— 

b . 92 

— 

b 2. 15 

55 

1.99 

.90 

1. 09 

3. 16 

2. 34 

108 

2. 30 

.81 

1. 5 

4. 16 

2.46 


Coefficient in eq. (14). 
^Extrapolated value. 


TABLE HI. - COMPUTED RESULTS BASED 


ON LANGMUIR PROBE DATA 


Magnetic field, 
B, 

tesla 

a l’ a 

V/m 2 

•t a 

b r 

V 

Electric field 
at r = a, 
V/m 

28X10" 4 

960 

33. 95 

96 

55 

2240 

42. 3 

224 

108 

5120 

43. 1 

512 


Coefficient in eq. (17). 
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0 20 40 60 80 100 120xl0" 4 

Magnetic field, B, tesla 

Figure 2. - Floating potential as a function of magnetic field. 



Figure 3. - Plasma noise for magnetic field value of 50xl0" 4 tesla and a floating 
anode; x = 2 microseconds per centimeter, y = 10 millivolts per centimeter. 
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5.00 3.75 2.50 1.25 0 1.25 2.50 2.75 5.00 


Position, r, cm 

Figure 4. - Molybdenum button front profile for various values of magnetic field. Retrace 
from X-Y plotter. 
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5.00 3.75 2.50 1.25 0 


Position, r, cm 

Figure 5. - Molybdenum button back profile for various values of 
magnetic field. Retrace from X-Y plotter. 



Position, r, cm 


Figure 6. - Molybdenum button data from Langmuir probe for vari- 
ous values of magnetic field. Retrace from X-Y plotter. 
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Impingement current, mA Number den 


Cf\ 

e 

to 


JA 1U 



1.6 

1.4 — 

1.2 


1.0 

.8 




.6 


0 10 20 30 40 30 60 70 80 90 100 110 120x10* 

Magnetic field, B, tesla 


Figure 7. - Anode impingement current and center plane number 
density as a function of magnetic field. 



Figure 8. - Langmuir probe trace. Magnetic field, 108xl0“ 4 tesla; 
radial position, 3.75 centimeters. Retrace from X-Y plotter. 


Current, arbitrary units 



Figure 9. - Log plot of Langmuir probe data 
from figure 8. Magnetic field, 108x10 " 4 
tesla; radial position, 3. 75 centimeters. 



Figure 10. - Plasma potential as a function of 
radial position for various magnetic fields. 
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Figure 11. - Calculated classical electron mo- 
bility as a function of magnetic field (eq. (5)) 
for constants appropriate to present experi- 
ment. Neutral number density, 2x10^ 
meter"^; average cross section, 14xl(T* 8 
square meters; average velocity, 10 6 meters 
per second. 



Figure 12. - Comparison of experimental values to theo- 
retical Bohm's equation (with Bohm’s constant 
C B = 0. 12(m 2 HT)/(V)(sec». 
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Figure 13. - Correlation study for 88 kilohertz per second. Magnetic field, 

55xl0" 4 tesla; x - 5 microseconds per centimeter, y - 10 millivolts per 
centimeter; probe phase shift, 120°. 



Figure 14. - Correlation study for 125 kilohertz per second. Magnetic field, 

45xI0~ 4 tesla; x - 2 microseconds per centimeter, y = 10 millivolts per 
centimeter; phase shift, 125°. 
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Figure 15, - Correlation study for 140 kilohertz per second. Magnetic field, 
97xl0" 4 tesla; x » 2 microseconds per centimeter; y = 10 millivolts per 
centimeter; phase shift, 136°. 


Figure 16. - Correlation study for 140 kilohertz per second. Magnetic field, 

97x10 ~ 4 tesfa with magnetic field reversed; x » 2 microseconds per centi- 
meter; y = 10 millivolts per centimeter; phase shift, 117°. 
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Magnetic field, B, tesla 

Figure 17. - Rotational frequency as a function of mag- 
netic field strength. 
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